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APOLLO MISSION D 
PERFORMANCE ANALYSIS 
OF RENDEZVOUS CHARTS 




1.0 Summary 

A digital analysis of the D Mission backup rendezvous 
charts has been performed to verify their ability to pre- 
dict CSI, CDH, TPI, and midcourse corrections in the presence 
of system and application errors and trajectory dispersions. 
Procedures for data acquisition were as defined in the cur- 
rent LM Rendezvous Procedures Document (Reference l). The 
charts simulated were those from the LM 3 onboard data 
package. The mission situation simulated represented the 
PNGS inoperative, rendezvous radar information available 
from the tapemeter, attitude data from the AGS and control by 
the CES. However, use of the CSI and CDH charts is presently 
ground ruled out by Reference 2. The study showed that the 
standard deviation of the arrival time at TPI was 5.5 minutes. 

The mean AT from TPI to intercept with braking and line-of-sight 
control was 35*8 minutes with a standard deviation of 1.5 minutes. 
If no braking or line-of-sight control was executed after the 
second midcourse correction, the mean miss distance would 
have been 0.6 n.m. with a standard deviation of .35 n.m. 

The mean total translational AV required after insertion was 
158 ft/sec with a standard deviation of 20 ft/sec. 
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2.0 Discussion 

This analysis was performed to determine the ability of 
the charts to predict all maneuvers after insertion required to 
complete the LM-active D mission rendezvous. One hundred inde- 
pendent sets of initial conditions (iC’s) were selected by 
adding random errors to the 6 components of the nominal relative 
state vector between the LM and CSM. Half of the, relative 
error was incorporated into the inertial state vector of each 
vehicle. 

The 100 cases were run including system and execution 
errors with braking and 13 ne -of- sight (LOS) control. The 
runs were repeated omitting braking and LOS control to obtain 
miss distance. The results of the 100 Monte Carlo runs were 
processed to obtain statistical data for the parameters of 
interest. 

The f3rst 25 cases were also run without system or 
application errors to establish the theoretical chart capabilities. 
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2.1 Study Rules 

The following ground rules consistent with D Mission 
procedures and planning were used in this study: 

1. CSX occurs 40 min 5 6 sec after insertion. 

2. CDH occurs 44 min 2 6 sec after C8I* 

3. TPI was assumed to occur 8 minutes after the 
elevation angle of 19*73 degrees was reached. 

4. The first midcourse correction (MCI) occurred 
10 minutes after TPI and the second midcourse 
correction (MC2) occurred at 22 minutes after 
TPI. 

5. CSI and CDH were burned in local vertical 
coordinates using impulsive thrust. 

6. TPI, MCI, MC2, braking, and LOS control 
were executed along and normal to the line- 
of-sight using finite thrust, burning each 
component individually. 

7. 130 degrees of CSM orbit travel between TPI 
and TPF. 

8. No out-of -plane corrections were made prior to 
LOS control during the braking phase. 
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2.2 Digital Program 

In the analysis, functions describing the backup 
rendezvous charts were programmed into a digital routine 
which integrated the equations of motion of two particles 
about an oblate planet. These functions allowed simulation 
of the procedures for using the backup rendezvous charts by 
incorporating elevation angle (AGS address 30k) 9 range, 
and range rate from the tapcmeter at the times called for 
by the backup data sequence. System errors were included 
in the data taken at each point and an appropriate error 
of application was included in each maneuver. 

The runs included effects of both bias and random 
errors as defined in Section 2.5. Bias errors were selected 
by the program at the beginning of each run and held constant 
for that run. Random errors were selected at each point 


data were taken. 
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2.3 Charts 

The charts modeled in the study were those which will be 
carried on the D Mission, and are shown in Figures 1-6. The data 
sequence was obtained from Reference 1 and is summarized on the 
relative motion plot of the nominal trajectory from CSX-40 to TPF 
(Figure 7). 

2.3.1 CSI 

The CSI chart solution is based on a Maclaurin’s expansion 
of four variables for the delta V at CSI (range rate at 30, 20, 
and 10 minutes prior to CSI and range 10 minutes prior to CSI). 

The coefficients are determined by the simultaneous solution of 
several expansions, each representing a dispersed trajectory prior 
to CSI. 

2.3.2 CDH 

The CDH chart solution utilizes the sinusoidal time history 
of range rate variations from coellipticity and relative velocity 
errors from coellipticity. Range rate data for the CDH chart are 
taken 29, 18, and 7 minutes prior to CDH. 

2.3.3 TP I and MC C 

The TPI charts solve for the relative position and velocity at 
TPI resolved into normal and along the line-of~sight (LOS) coordinates. 
The measured relative conditons are differenced from the required 
conditions for intercept in 130 degrees of orbit travel. Informa- 
tion required for the TPI charts is the elevation angle at 8 and 
5 minutes prior to TPI and range and range rate at 5 minutes before 
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TPI. In a similar manner, the MCC charts maintain the time of 
TPF consistent with the TPI maneuver. Data are taken for the 
midcourses 5 and 8 minutes after TPI for MCI and at 17 and 20 min 
utes after TPI for MC2. The same measurement sequence as used 
at TPI is used for both midcourse corrections. 
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» 

2.4 Initial Conditons 

. The initial conditions were generated by perturbing the 
nominal state vectors of the LM and CSM with errors supplied by 
a relative covariance matrix. Half the relative error was applied 
to the state vector of each vehicle. 

The nominal vectors relative to the LM orbital plane were 
derived from Reference 3 and are summarized as follows: 

Table 2-1 
Initial Conditions 
Insertion + 56 Sec 
(95 hrs 4l min 48 sec) 


LM Altitude 

863730 feet 

CSM Altitude 

798403 feet 

LM Total Velocity 

25431.2 ft/sec 

CSM Total Velocity 

25468.7 ft/sec 

LM Flight Path Angle 

-.00563 deg 

CSM Flight Path Angle 

-.00652 deg 

In Plane Central Angle 

.401 deg 

LM Latitude 

1.4 deg 

Out of Plane Distance 

-31.3 feet 

Out of Plane Velocity 

-.02 ft/sec 

Heading Angle (relative to equator) 

29 deg 



The covariance matrix used for initialization of 


the Monte Carlo runs is as follows: 

Table 2-2 
Covariance Matrix 


9725775. -358310.0 

11.8853 

799.127 -205.801 - 

-.043464 

-358310.0 12124470. 

-6.36828 

-887.559 199.080 

.0249741 

11.8853 

-6.36828 

6262300. 

.0189843 -.004047 

20.2459 

799.127 

-887.559 

.0189843 

11.35113 -.42248 • 

-.0000708 

-205.801 

199. 080 

-.004047 

-.42248 3.202925 

.00001599 

-0434646 

.029741 

20.2459 

-.0000708 .00001599 2.394625 


It was obtained by increasing between 9 and 25 times the 
diagonal elements of a post insertion covariance matrix 
provided by Math Physics Branch. 
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2.5 Errors in Sensors and Execution 

The la errors in sensor and maneuver executions were; 

1. System Errors 

A. Noise 

1 ) Range . 333# 

2) Range Rate .433# or .433 ft /sec 

which ever is larger 

3) Elevation Angle .12 degree 

B. Biases and Drifts (constant for a given run) 

1) Range Rate .333 ft /sec 

2) Initial Pitch Bias .1 deg (assumes 

calibrated COAS) 

3) Pitch Drift Rate .23 deg/hr 

2. Execution Errors 

A. Reading Tapemeter 

1) Range Rate .25 ft /sec 

2 ) Range 

a) Outer Scale 2400 ft 

b) Middle Scale 100 ft 

B. Application of Burns .5 ft /sec (per axis) 


C . Time Measurements 


1.0 sec 
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2,6 Braking Schedule and LOS Angular Rate Deadbands 

The braking schedule used in this simulation consisted 
of five gates arid a lower limit on the range rate. The first 
gate was at 13500 feet. At this point only LOS control was 
executed because the allowed range rate was 80 ft /sec. The 
second gate was at 6000 ft with an allowed range rate of 
30 ft /sec. The nominal range rate at this range was 29 ft /sec. 
The remaining gates were 20 ft/sec at 3000 ft, 10 ft/sec at 
1500 ft and 5 ft /sec at 500 ft. The lower range rate limit 
consisted of a straight line connecting 20 ft /sec at 13500 ft 
and 0 ft /sec at intercept. 

Both the upper and lower range rate limits are shown in 
Figure 8. LOS control procedures were simulated by sampling 
inertial drift of the LOS inplane and normal to the orbit 
plane every 15 seconds beginning at a range of 13500 ft. 

When the LOS rates exceeded .3 mr/sec at a sampling time, 
thrust was applied in the appropriate axis in increments 
of 1 second until the LOS rate was reduced below the 
threshold. The 15 seconds were allowed to elapse before 


sampling again. 
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3.0 Results 

Several sets of Monte Carlo runs were made to obtain 
statistical data for determination of the Effects of errors, 
trajectory dispersions, and braking on the size of maneuvers, 
arrival time at TPI, and total translational AV. The sets 
of runs are identified in the following table: 




Table 3-1 
Run Summary 



Set 

Number 
of Runs 

Errors in 

Maneuver 

Solutions 

Maneuvers 

Applied 

LOS Control 
and Braking 

A 

100 

Yes 

With Errors 

Yes 



No 



B 

100 

Yes 

With Errors 

No 

C 

LTV 

CVi 

No 

No Errors 

No 

D 

25 

No 

No Errors 

Yes 


Solutions 

for the maneuvers 

in SET A were 

obtained with 


and without errors so that the effect of sensor and reading 
errors on the chart solutions could be determined. However, 
all maneuvers for SET A were made usi..: r ; the solutions with 
errors. The runs for SET B were identical to SET A, but 
with braking and LOS control ommitted to establish 
miss distance. 

SETS, C and D were run to establish baseline data for 
chart performance. It was felt that a reduced number of runs 
would suffice to obtain statistically meaningful results since 
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only initial conditions were varied. Examination of sig- 
nificant parameters such as maneuver solution and miss 
distances revealed nearly normal statistical distributions, 
confirming the adequacy of the 25 runs for those sets. The 
runs for SETS C and D were made with the same initial con- 
ditions as the first 25 runs of SET A’. 
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3»1 Maneuver Value 3 

The nominal chart solutions, average, mean, and standard 
deviation for each maneuver in SET A and SET D are shown in 
Table 3-2 on Page 3-^. The data given for SET A are the 
solutions with errors. The average, mean, and standard devia- 
tion for the difference between the error solution and the no 
error solution computed for each maneuver in SET A are shown 
in Table 3-3 on Page 3-5 • 

It can be noted from the data on Tables 3-2 and 3-3 that 
the chart solutions with errors progressively decrease in 
accuracy from CSI to CDH to TPI. The trend then reverses with 
MCC1 more accurate than TPI and MCC2 being the most accurate 


of the maneuvers 
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Table 3-3 

Differences Between Chart Solutions 
With Errors and Without Errors in SET A 


Maneuver 

Average 
ft/ sec 

Mb an 
ft/seo 

Standard Deviation 
ft /sec 

CSI 

.79 

.01 

1.02 

CDh^dVy 

• 94 

*07 

1.24 

CDH ^V H 

.48 

• . 06 

*59 

tpi ^ v los 

2.31 

-.04 

2.89 

TFI <UV N 

2.06 

-.08 

2.57 

MCCi4av L0 s 

2.50 

-d- 

CVl 

• 

1 

3.07 

MCCl/i4V N 

1.19 

-.14 

1.45 


•47 

-.00 

*59 

MCC2A1 V n 

.54 

.02 

.66 


The data of Table 3-2 indicates the amount of error 
directly attributable to the sensor and reading errors 
listed in Section 2.5* 
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3.2 Mias Distance 

The miss distance was established by sets B and C. The 
average in plane miss distance at the point of closest 
approach for the 25 cases without errors (SET C) was 
513 feet, and for the 100 cases with errors (SET B) was 
2505 feet. 

The average, mean, and standard deviation of the 
components of the miss distance in a local vertical 
coordinate system with X along the radius vector of the 
LM, Z along the angular momentum vector of the LM, and 
Y completing the right-handed system were as follows: 

Table 3-4 

Coordinates at Closest Approach 


Axis Average Mean Standard Deviation 



SET B 

Feet 

SET C 

Feet 

SET B 

Feet 

SET C 
Feet 

SET B 

Feet 

SET C 
Feet 

X 

1482 

375 

-915 

-375 

1802 

219 

Y 

2020 

350 

-802 

-350 

2710 

188 

Z 

1844 

1508 

-12k 

- 95 

2332 

1879 
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3.3 AV Used 

The mean total translation AV used in the 25 cases without 
errors (SET D) was 138.2 ft/sec with a standard deviation of 
13.1* ft/sec while the mean for the 100 cases with errors (SET A) 
was 158.7 ft/sec with a standard deviation of 20.5 ft/sec. 

The minimum und maximum AV oases without errors required 
108.2 ft/sec and 170.8 ft/sec, respectively, while with errors 
minimum and maximum AV cases required 100.0 ft/sec and 227.0 ft/sec. 
Figure 9 shows the distribution of total AV. 

A breakdown of how the AV was used is shown in the following 
table : 


Table 3-5 
Translation AV 


Maneuver AV (SET D) 

Average Without Errors 
ft/sec 


AV (SET A) 
Average With Errors 
ft/sec 


CSX 

37.1* 

36.8 

CDH 

40.9 

1 * 0.2 

TFI 

22.1 

2 l *.0 

MCC1 

1.9 

6.8 

MCC2 

1.3 

10.8 

Braking and LOS 
Control 

35.2 

1 * 0.0 
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3.U Arrival Time at TP I 

The mean arrival at TPI for the 25 cases without errors 
(SET C) was 5 seconds later than nominal with a standard de- 
viation of 18 seconds, while the mean fop the 100 cases with 

errors (SET A) was 12 seconds late with a standard deviation 
of 5 min and 35 seconds. 

Figure 10 shows the distribution of arrival’ time at TPI 
over intervals of two minutes for the 100 cases with errors. 
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3.5 A T from TPI to TPF 

The mean AT from TPI to close approach without braking for 
the 25 cases in SET C was 1980 seconds with a standard deviation 
of 18 seconds while the mean AT of transfer without braking for 
the 100 cases in SET B was 1963 seconds with a standard devia- 
tion of 72 seconds. 

The mean AT from TPI to intercept with braking and LOS 
control for the 25 cases in SET D was 2121 seconds with a standard 
deviation of 22 seconds while the mean AT of transfer with 
braking for the 100 in SET A was 2148 seconds with a standard 
deviation of 87 seconds. The nominal case required 2115 seconds 
with braking and LOS control. 

Figure 11 shows the distribution of the AT transfer with 
braking for the 100 cases with errors. 
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TERMINAL PHASE ^ INITIATION 




LW-120 CSM-130 

130 degree transfer FIRST MIDCOURSE CORRECTION 

9-17 68 7 JET Z * JET X 
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